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Nomenclature 


a  =  ‘quality’  indicator  of  ionization  in  the  ®Vdf  contribution  of  the  O-Factor,  a  =0.99 

Ej  =  first  experimental  parameter  group  =  !4(T /  m  )2/(V^/M)  =  <D(l-(3)x. 

E2  =  second  experimental  parameter  group  =  ( lit  /Ia)(5/J4 )  =  ( l-r)/'/_. 

5  =  Faraday  constant,  96485  coulombs/mol  of  charge. 

fj  =  ion  mass  fraction  at  exit=  fi+f2+f3- 

f*i,f* 2,f* 3  =  reduced  ion  mass  fractions  at  exit,  f* i+f* 2+f* 3=1. 

g0  =  Earth’s  gravitational  constant  at  sea  level,  9.806  m/s2. 

Ia  =  measured  anode  current. 

Isp  =  specific  impulse,  independent  of  unit  system,  Isp  =  <v>/gG. 

lit  =  measured  propellant  anode  mass  flow  rate, 

m  =  measured  propellant  cathode  mass  flow  rate. 

M  =  molecular  weight  of  propellant,  xenon  =  0. 13 13  kg/mol. 

Pin  =  input  power  to  anode. 

Pjet  =  power  of  jet  =  ’A  m  <v2>. 

P magnets  =  input  power  electromagnets. 

Q  =  average  charge  of  xenon  ions,  f  i  +  2f*2  +  3f*3. 

q  =  unit  of  charge,  1.609  x  10"19  C. 

r  =  fractional  loss  of  current,  electron  recycle  fraction. 

(1-r)  =  current  utilization  efficiency. 

T  =  measured  thrust  vector. 

Va  =  measured  anode  voltage. 

<v>  =  mass  weighted  average  exit  velocity,  specific  impulse. 

<v2>  =  mass  weighted  average  squared  exit  velocity. 

<  v2>,  =  mass  weighted  average  squared  ion  exit  velocity, 

a  —  plume  momentum  divergence  half-angle. 

P  =  fractional  loss  of  acceleration  voltage. 

( 1  -(3)  =  voltage  utilization  efficiency. 

AV  =  average  acceleration  voltage  =  (f* i  5Vi+  2f* 2  5Vi+  3f* 3  6Vj)  /  Q. 

5V!,  5V2,  5V3  =  acceleration  voltages  of  Xe+1,  Xe+2,  Xe+3. 

P energy  =  energy  efficiency  =  Pjet/Pi„. 

r| thrust  =  experimental  thrust  efficiency. 

h total  =  total  thruster  efficiency  =  ii,hrust  [  lit  /( lit  +  mc )]  [Pin  /(  Pin  +  Pmagnets)] 

(f)  =  O-Factor  =  Ovdf  ‘Pdiv  =  <v>2/<v2>,  propellant  utilization  efficiency. 

Ovdf  =  velocity  distribution  component  of  the  O-Factor,  =  a  fj 

X  =  X_Factor  =  fQ,  average  charge  per  particle  (ions  and  neutrals),  beam  current  per  mol  of  propellant. 


I.  Introduction 

HALL  thrusters  utilize  orthogonal  electric  and  magnetic  fields  to  electrostatically  accelerate  ions  through  a 
discharge  channel,  while  electrons  travel  in  an  azimuthal  trajectory  around  the  channel  annulus  with  a  small 
axial  velocity  component  leading  to  the  anode.  Complex  coupling  of  plasma  properties  and  perturbation  effects 
make  physical  parameters  such  as  beam  divergence,  ionization  fraction,  and  recycled  electron  current  difficult  to 
measure  directly.  Flowever,  these  values  can  be  extracted  using  macroscopic  performance  and  far-field  plume 
measurements.  An  investigation  of  thruster  operation  at  constant  anode  flow  (20  mg/s)  and  constant  power  (6  kW) 
is  examined  to  demonstrate  the  utility  of  thruster  utilization  efficiencies  and  the  efficacy  in  bounding  a  consistent  set 
of  plasma  properties  from  a  limited  number  of  macroscopic  measurements. 
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II.  Hall  Thruster  Efficiency  Analysis 

A  Hall  thruster  efficiency  architecture  developed  at  AFRL  analytically  separates  thruster  efficiency  into  the 
product  of  three  utilization  efficiencies:  voltage,  current,  and  propellant  utilization  efficiency.  This  method 
accounts  for  the  effects  of  the  ion  velocity  distribution  function  (VDF),  ionization  fraction,  multiply  charged  ion 
species,  electron  recycle  fraction,  plume  divergence,  and  the  loss  of  acceleration  voltage  within  the  discharge 
chamber.  The  architecture  is  formulated  such  that  energy  efficiency  becomes  naturally  expressed  as  a  product  of 
voltage  and  current  utilization  efficiencies,  and  is  rigorously  separated  from  propellant  utilization  efficiency.  Thus, 
thrust  efficiency  is  the  product  of  propellant  utilization  efficiency  and  energy  efficiency.  A  detailed  account  of  the 
following  efficiency  architecture  with  a  comprehensive  examination  of  Hall  thruster  physics  and  its  manifestation  in 
thruster  utilization  efficiencies  is  detailed  in  the  work  of  Larson,  Brown,  and  Hargus. 1 

Hall  thruster  efficiency  is  defined  as  the  conversion  of  anode  electrical  energy  to  directed  momentum  to  produce 
thrust,  and  experimentally  determined  as: 
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thrust 


1  •  2 
m  <  v  > 

2 


1  {  T 


21  P, 


<  v  > 


(1) 


In  Eq.  (1),  thrust  efficiency  may  also  be  written  as  the  product  of  the  thrust-to-power  ratio  (T/Pin)  and  <v>,  where 
<v>  is  the  propellant  jet  specific  impulse  (impulse  per  unit  mass)  with  dimensions  of  velocity  in  SI  units.  This 
expression  reveals  the  trade-off  between  T/P;n  and  specific  impulse  for  a  given  thruster  efficiency,  and  explicitly 
shows  that  increasing  specific  impulse  to  high  values  at  constant  efficiency  requires  diminishing  thrust  and/or 
increasing  power.  Further  analysis  of  the  thrust  efficiency  provides  insight  into  losses  in  Hall  thruster  performance. 


thrust 


<V> 


? 

<  V-  > 


m  <  v'  > 


I  V 

a  a 


=  0-^  =  0  q 

p  'energy 


where 


(2) 


q>  = 


<  v>" 


2 

<  V  > 


(3) 


The  derivation  in  Eq.  (2)  separates  thrust  efficiency  into  the  product  of  energy  efficiency  and  the  (D-Factor. 
Energy  efficiency  is  the  ratio  of  jet  kinetic  energy  to  anode  electrical  energy.  The  cD-Factor  defined  in  Eq.  (3) 
represents  propellant  utilization:  it  is  unity  for  100%  ionization  to  a  single  ion  species  where  all  thrust  vectors  are 
directed  along  the  same  axis.  It  is  a  direct  relation  between  the  ion  momentum  and  kinetic  energy,  which  accounts 
for  losses  due  to  incomplete  ionization,  plume  divergence  and  non-uniformity  of  the  VDF  that  occurs  when  the  jet  is 
composed  of  multiple  ion  species  with  varying  velocities. 

Energy  efficiency  can  be  further  divided  into  the  product  of  voltage  and  current  utilization.  The  y- factor  is 
introduced  in  Eq.  (4)  to  account  for  the  average  particle  charge  of  propellant  in  the  plasma,  including  unionized 
neutral  propellant.  It  is  defined  as  the  product  of  particle  ionization  fraction,  f.  and  average  ion  charge,  Q,  in  Eq.  (5) 
for  a  plasma  with  three  ion  species. 
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(5) 


Although  the  y-factor  cancels  in  the  energy  efficiency  term,  it  reveals  the  physical  nature  of  the  voltage  and 
current  utilization.  The  voltage  utilization,  (1-P),  shown  in  Eq.  (6)  compares  the  average  potential  energy  gained  by 
ions  with  the  total  potential  available  from  the  anode  discharge  voltage. 
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To  simplify  the  analysis,  all  ions  are  assumed  to  be  created  at  the  same  location.  Thus,  the  potential  available  for 
acceleration  is  AV  and  the  voltage  utilization  efficiency  is  (1-P)  =  AV/Va.  This  utilization  is  calculated  from  the 
most  probable  ion  voltage  with  an  RPA.  Voltage  utilization  is  maximized  for  complete  ionization  at  the  anode  face. 
The  approximation  used  in  Eq.  (6),  in  which  average  neutral  to  ion  specific  kinetic  energy  is  equal  to  f;  (i.e., 
<v2>/<v2>+  ~  fp,  is  accurate  to  better  than  99.999%.  This  calculation  is  based  on  the  velocity  ratio  of  singly  charged 
xenon  ions  to  unionized  neutral  propellant  being  greater  than  or  equally  to  100,  which  was  experimentally 
determined  using  electrostatic  probes  by  Elargus,  Meezan,  and  Cappelli  for  a  low  power  Hall  Thruster. 

The  total  ion  beam  current  exiting  the  thruster  is  given  by  [y  lit  (ffjM)].  The  ratio  of  ion  beam  current  to  the 
discharge  current  is  the  current  utilization  efficiency  (1-r),  as  defined  in  Eq.  (7).  The  quantity  ‘r’  is  the  electron 
recycle  fraction,  and  represents  the  electrons  emitted  from  the  cathode  that  are  not  used  for  beam  neutralization. 


(1  -r)  =  —  —  x  =  —  —  fj(fi*  +  2f>*  +  3f3* ) 
L  M  L  M  lW  3 


(7) 


Current  utilization  approaches  unity  as  the  number  of  recycled  electrons  approaches  zero.  Elowever,  there  is  a 
minimum  value  of  ‘r’  required  to  sustain  ionization  processes  in  the  thruster  and  provide  the  energy  necessary  for  a 
given  ionization  fraction  and  ion  species  distribution. 

Two  experimental  parameter  groups,  Ej  and  E2,  may  be  written  using  the  preceding  utilization  efficiencies. 
These  parameters  are  calculated  based  solely  on  thruster  telemetry  and  performance  measurements,  and  provide 
insight  about  the  relative  magnitudes  of  the  individual  utilization  efficiencies.  The  product  of  Ej  and  E2  is  equal  to 
thruster  efficiency  as  shown  in  Eq.  (8). 
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Equations  (9)  and  (10)  define  how  Ei  and  E2  depend  on  the  utilization  efficiencies.  The  quantity  Ej  is  a  function 
of  propellant  utilization,  voltage  utilization,  and  the  y-factor.  It  is  important  to  note  that  since  both  the  y-factor  and 
O- Factor  are  directly  proportional  to  the  ionization  fraction,  the  experimental  parameter  E,  is  proportional  to  f2. 
Since  the  quantity  E2  is  inversely  proportional  to  f.  it  is  expected  that  E|  will  be  more  sensitive  to  changes  in  the 
ionization  fraction. 
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In  the  absence  of  diagnostics  for  the  determination  of  plasma  properties,  these  experimental  parameter  groups 
allow  limits  to  be  placed  on  acceptable  values  for  the  average  charge  Q,  ionization  fraction  f .  and  momentum 


4 

The  30th  International  Electric  Propulsion  Conference,  Florence,  Italy 
September  17-20,  2007 


divergence  half-angle,  a.  Analysis  using  Ej  will  ignore  variations  in  the  quality  term  ‘a’  from  Eq.  (9),  since  this 
term  is  approximately  0.99  for  most  thrusters  and  rarely  deviates  more  than  2%  from  unity.  A  hill  explanation  is 
found  elsewhere  in  the  literature.13 


III.  Experimental  Apparatus 


A.  AFRL  Test  Facilities 

Performance  verification  tests  at  AFRL  were  conducted  in  Chamber  3,  a  stainless  steel,  cylindrical  vacuum 
chamber  3.3  m  diameter  by  8-m  long.  The  thruster  exit  plane  is  located  in  the  center  of  the  chamber  diameter, 
pointing  down  the  long  axis  approximately  7  m  from  the  far  end.  All  surfaces  that  are  optically  visible  to  the  plasma 
plume  are  shielded  with  14”  thick,  high  purity,  sulfur-free  carbon  plates.  A  dedicated  dry  Stokes  mechanical  pump 
reaches  rough  vacuum  with  a  pumping  speed  of  450  L/s  and  a  Pfeiffer  Vacuum  turbomolecular  pump  removes 
lighter  gases.  The  facility  utilizes  eight  helium-cooled  cryopanels,  resulting  in  a  maximum  xenon  pumping  speed  of 
140,000  L/s.  The  heat  load  to  these  cryopanels  is  moderated  by  shrouds  coated  with  low  emissivity  paint,  which  are 
cooled  by  four  IGC  Polycold  Cryogenic  Refrigeration  units.  A  conservative  measurement  of  the  ultimate  base 
pressure  is  less  than  1  x  10'7  torr. 

An  MKS  Instruments  cold  cathode  ionization  gauge,  located  on  the  chamber  ceiling  above  the  thruster, 
registered  background  pressure.  Using  a  Xenon  calibration  factor  of  2.87,  the  adjusted  pressure  for  flow  rates  of  10, 
15,  and  20  mg/s  were  8  x  10"6,  1.3  x  10’5,  1.5  x  10'5  Torr.  High  background  pressure  artificially  raises  performance 
due  to  ingestion  of  background  neutrals  into  the  discharge  chamber.  Randolph’s  criterion  asserts  that  facility  effects 
have  negligible  influence  on  performance  for  background  pressures  below  2x  1 0"  Torr.2  Since  all  performance 
measurements  were  taken  at  pressures  below  this  threshold,  it  was  determined  that  neutral  entrainment  at  these 
pressure  and  flow  rates  had  a  negligible  impact  on  thruster  performance.3  4  No  corrections  to  performance 
measurements  were  made  to  compensate  for  facility  effects. 

The  propellant  flow  system  was  supplied  with  research  grade  Xenon  (99.999%  purity)  via  %  inch  stainless  steel 
tubing.  Anode  flow  was  controlled  with  a  1000  seem  Unit  1661  digital  mass  flow  controller.  Cathode  flow  was 
regulated  with  a  50  seem  Unit  7300  digital  mass  flow  controller.  All  mass  flow  controllers  were  calibrated  with 
Xenon  prior  to  testing  with  a  DH  Instruments  Molbloc/Molbox  gas  flow  calibration  system  with  an  estimated 
uncertainty  of  ±0.5%  according  the  manufacture’s  specifications. 

Thruster  discharge  power  was  provided  by  two  Electronic  Measurement  Inc.  ESS  Power  Supplies  arranged  in 
parallel,  each  capable  of  600V,  16A  output.  A  simple  RLC  circuit  provided  low  pass  filtering  of  the  main  discharge, 
and  the  thruster  body  was  connected  to  chamber  ground.  Discharge  voltage,  discharge  current,  and  cathode 
potential  were  monitored  through  high  accuracy  current  shunts  with  an  uncertainty  of  ±0.5%.  Coils  were  energized 
with  Sorenson  DLM  60-10  power  supplies.  The  cathode  heater  was  powered  with  an  EMS  20-30  supply  and  the 
cathode  keeper  used  a  Sorenson  LHP  300-3.5  power  supply.  The  heater  power  supply  was  shut  down  during 
thruster  operation. 

Thrust  measurements  were  taken  with  a  null  type,  inverted  pendulum  thrust  stand.  At  steady  state  conditions, 
thrust  measurements  varied  by  <0.5  mN.  A  diluted  Glycol  solution  at  10°C  cooled  the  copper  shroud  surrounding 
the  thrust  stand  to  decrease  thermal  drift.  Thrust  stand  calibrations  were  performed  in  several  configurations  to  test 
the  influence  of  electromagnet  coils  and  cathode  operation.  Repeated  tests  showed  consistent  thruster  performance 
to  within  1%,  signifying  that  the  calibrations  and  zero  offsets  were  unaffected  by  these  variables.  As  the  thrust  stand 
and  thruster  reached  thermal  equilibrium,  the  calibration  curves  showed  repeatability  with  less  than  1%  variation. 
The  zero  offset  was  recorded  after  each  test  point  and  used  with  the  calibration  slope  to  calculate  thrust.  In  cases 
where  the  thrust  stand  was  not  calibrated  after  each  operating  condition,  the  zero  position  was  found  for  each  point 
from  a  linear  interpolation  of  the  zero  offset  before  firing  and  after  shutdown.  Inclination  was  managed  with  a 
stepper  motor  to  level  the  thruster  and  minimize  null  coil  current  during  high  thrust  operation.  This  minimized 
thermal  drift  in  the  null  coil  signal  after  shutdown  and  yielded  a  more  reliable  zero  offset.  Due  to  the  highly 
repeatable  thrust  measurements  and  minimal  thermal  drift,  uncertainty  in  thrust  is  estimated  at  ±1%. 

A  22-bit  Agilent  Datalogger  transmitted  thrust,  mass  flow,  and  power  to  a  LabView  code  monitoring  the  anode 
efficiency  and  thruster  performance  in  real-time.  For  a  given  discharge  voltage  and  flow  rate,  the  coil  currents  were 
adjusted  for  maximum  anode  thruster  efficiency.  When  peak  efficiency  was  achieved,  the  magnet  coils  were  further 
tuned  for  maximum  T/P  ratio.  The  efficiency  optimization  did  not  always  correspond  to  minimum  discharge 
current,  since  decreasing  thrust  would  further  reduce  anode  thruster  efficiency  (Efficiency  ~  Thrust2  /  Current). 
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B.  Hall  Thruster 

A  6  kW  laboratory  model  test-bed  Hall  thruster 
was  used  in  this  investigation.  A  large  thruster  was 
utilized  to  enable  placement  of  a  centrally  mounted 
internal  cathode  that  was  directed  along  the 
thruster  centerline.  The  schematic  in  Fig.  1 
illustrates  the  internal  cathode  configuration  in  a 
typical  Hall  thruster  design.  The  thruster  utilizes 
inner  and  outer  electromagnets  to  form  standard 
magnetic  lens  topography.8  Constant  anode 
propellant  flow  operating  conditions  at  20  mg/s 
ranged  from  100  -  400  V  discharge.  Constant 
power  operating  conditions  at  6  kW  ranged  from 
200  -400  V.  Cathode  flow  fraction  was  7%  at  all 
operating  points. 


Recycled  e 


Beam  Neutralizing  e 

THRUSTER  CENTERLINE 


C.  Plasma  Diagnostics 

Plume  examinations  were  executed  one  meter 
downstream  of  the  thruster  exit  plane  as  a  function 
of  angular  position  from  centerline.  The  plasma 
diagnostic  array  at  AFRL,  shown  in  Fig.  2,  includes 
a  retarding  potential  analyzer  (RPA),  Langmuir 
probe,  gridded  Faraday  probe,  and  nude  Faraday 
probe.  A  rotational  arm  with  axis  of  rotation  located 
above  the  exit  plane  on  thruster  centerline  is 
controlled  with  a  Thermionics  Stepping  Motor 
Controller.  The  probe  array,  as  shown  in  Figure  3, 
is  swept  continuously  180  degrees  in  a  constant  one 
meter  radius  arc  for  Faraday  probe  measurements. 
The  continuous  rotational  motion  at  low  speeds 
minimizes  mechanical  arm  vibration  and  prolongs 
probe  life  by  decreasing  time  in  the  plume.  RPA 
and  Langmuir  probe  data  collection  is  taken  on 
thruster  centerline  one  meter  downstream  of  the  exit 
plane,  with  the  Langmuir  probe  3  cm  off  axis. 


Figure  1.  Schematic  of  a  centrally  mounted  internal 
cathode  configuration  in  a  typical  Hall  thruster  design. 


/  \ 


W(ru-  ' 


Figure  2.  Plasma  diagnostic  array  at  AFRL.  Probes  from 
left  to  right  include:  Langmuir  probe,  RPA,  nude  Faraday 
probe,  and  gridded  Faraday  probe. 


Plasma  Diagnostic  Array 


0°  \  ~  Thruster  Centerline 
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1.  Retarding  Potential  Analyzer 

An  RPA  measured  the  ion  voltage  distribution 
on  the  thruster  centerline  with  a  series  of  positively 
and  negatively  biased  grids.  The  positive  potential 
of  an  ion  retarding  grid  acts  as  a  high  pass  filter,  and 
only  permits  ions  above  the  grid  threshold  voltage 
to  reach  the  collector.  The  resulting  ion  voltage 
distribution  is  found  from  the  derivative  of  the 
current-voltage  relationship. 5 

A  series  of  grids,  as  shown  in  Figure  4,  filters 
ambient  electrons  and  selected  ions  before  they 
reach  the  copper  collector.  Grids  are  composed  of  a 
0.005”  thick,  photochemically  etched  316  SS  with 

an  open  area  fraction  of  approximately  38%  and  Figure  3  p,asma  diagnostic  array  on  the  rotational  arm  at 
separated  by  MACOR  spacers.  A  phenolic  sleeve 
inside  the  stainless  steel  body  electrically  isolates 

the  grids.  A  ceramic  cover  with  a  1  cm  diameter  opening  reduced  the  ion  flux  at  high  flow  rates  and  avoided  current 
saturation  of  the  collector  electronics.  The  outer  grid  exposed  to  the  plume  is  floating  to  minimize  the  perturbation 
between  the  probe  and  ambient  plasma.  Electron  repelling  grids  2  and  4  are  biased  to  -15  V  below  thruster  ground 
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by  a  BK  Precision  1760  Triple  Output  DC  power 
supply.  Grid  2  repels  electrons  from  the  plume 
while  grid  4  suppresses  secondary  electrons 
resulting  from  ion  impacts  with  the  collector.  The 
ion  retarding  grid  is  connected  to  a  Keithley  2410 
1100  V  Sourcemeter  and  swept  from  0  to  -500  V 
with  respect  to  ground.  A  picoammeter  reads  the 
collection  current,  which  is  recorded  as  a  function 
of  voltage  applied  to  the  ion  repelling  grid.  Further 
RPA  construction  and  operation  details  may  found 
elsewhere.6,7 

2.  Langmuir  Probe 

A  single,  cylindrical  Langmuir  probe  parallel  to 
the  centerline  plume  was  used  to  measure  plasma 
potential.  The  probe  consisted  of  a  0.2  mm 
diameter,  15  mm  long  Tungsten  wire  connected  to 
a  Keithley  2400  Sourcemeter.  Langmuir  probe 
data  were  collected  during  the  RPA  sweep  and 
analyzed  to  determine  the  plasma  potential  for 
RPA  voltage  distribution  calculations.  The 
Langmuir  probe  was  positioned  3  cm  and 
approximately  1  degree  off  the  centerline  axis. 
Prior  studies  concluded  the  difference  in  plasma 
potential  varies  by  only  a  few  volts  across  large 
angles.8  Due  to  the  proximity  of  the  Langmuir 
probe  to  the  RPA,  the  difference  in  plasma 
potential  between  the  Langmuir  probe  location  and 
the  RPA  on  the  thruster  centerline  was  estimated  at 
less  than  one  volt. 

A  potential  diagram  outlining  the 
correlation  between  measured  values  with  the  most 
probable  ion  voltage  and  loss  voltage  is  shown  in 
Fig.  5.  The  average  acceleration  potential,  AV,  is 
found  in  Eq.  (11)  as  the  difference  between  the 
most-probable  ion  voltage  from  RPA  traces  (Vrpa) 
and  plasma  potential  (VPiasma).  The  loss  voltage 
(Vloss)  in  Eq.  (12)  is  subsequently  calculated  from 
the  difference  between  discharge  voltage  (Va)  and 
average  acceleration  potential.  Voltage  utilization 
efficiency  can  then  be  determined  from  the  average 


Grid  1 :  Floating 


Figure  4.  Cross-sectional  schematic  of  the  RPA  showing 
ion  and  electron  repelling  grids. 
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Figure  5.  Potential  diagram  illustrating  the  relationship 
between  measured  voltages  (V„,  VC.G),  plume  measurements 
(Vrpa,  VPlasma)  and  calculated  potentials  (AV,  VLoss). 

acceleration  voltage  as  shown  in  Eq.  (13). 
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Bias  Voltage 
Power  Supply 


3.  Nude  Faraday  Probe 

The  nude  Faraday  probe  measures  directed  ion 
flux  to  a  19  mm  collector  electrode  bordered  by  a 
43  mm  diameter  annular  guard  ring.  This  guard 
ring  is  biased  to  the  same  negative  potential  as  the 
collector  and  forms  a  uniform,  planar  sheath  in 
front  of  the  collection  area.  A  simplified  circuit 
diagram  is  shown  in  Fig.  6.  The  bias  voltage  is  set 
to  -15  V  with  respect  to  ground  by  a  Sorenson  LT 
30-3  DC  power  supply.  Collected  current  is  found 
from  the  voltage  drop  across  a  40.4  ohm  shunt 
resistor,  and  corresponding  current  density  is 
calculated  from  current  divided  by  the  collector 
electrode  area. 

The  95%  charge  divergence  half-angle  (a95) 
was  calculated  separately  for  left  and  right  sides  of  the  plume, 


Guard  Ring 


Measured 

Shunt 

Voltage 


40.4  Cl 


o 


Figure  6.  Simplified  circuit  diagram  of  the  nude  Faraday 
probe. 


and  determined  the  location  where  the  signal 
integrated  to  95%  of  the  total  collected  voltage  from  centerline  to  ±90  degrees,  as  shown  in  Eq.  (14). 


±90 


±a9c 


0.95 


|  j(0)sln(0)d0  =  j"  j(0)sin(0)d0 


(14) 


Beam  current  is  calculated  by  integrating  current  density  in  the  hemisphere  about  thruster  centerline  as  shown  in 
Eq.  (15).  The  ratio  of  this  calculated  ion  beam  current  to  the  discharge  current  provides  the  current  utilization 
efficiency,  (1-r). 


+90 

=m2  \mmeye 
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Previous  studies  have  found  a  strong 
dependence  on  facility  backpressure,  resulting 
in  an  artificial  increase  in  current  density  at 
half-angles  approximately  30  degrees  beyond 
thruster  centerline.  This  dependence  is 
attributed  to  the  presence  of  charge  exchange 
(CEX)  ions,  which  account  for  a  significant 
fraction  of  ions  on  outer  edges  of  the  plume. 
9,10.11.12  Current  collected  on  the  outer  wings 
is  erroneously  high  and  results  in 
overestimation  of  the  total  beam  current  and 
divergence  angles.  Fig.  7  illustrates  increased 
current  density  on  the  wings  for  a  one  meter 
radial  sweep.  Results  have  been  normalized  to 
the  maximum  signal  value  for  comparison 
between  the  nude  and  gridded  Faraday  probes. 
A  corrected  nude  Faraday  probe  current  density 
is  found  from  the  linear  extrapolation  on  a  log 
scale  for  angles  greater  than  30  degrees  off 
centerline. 13  Although  results  from  the  gridded 
Faraday  probe  also  show  evidence  of  increased 
current  density  on  the  wings,  the  correction  is 
unnecessary  since  low  energy  CEX  ions  are 


Figure  7.  Normalized  current  density  from  nude  and 
gridded  Faraday  probe  traces  showing  increased  current 
collection  on  the  wings  due  to  CEX  ions.  The  dashed  line 
illustrates  the  linear  extrapolation  of  nude  Faraday  probe 
data  on  a  log  scale. 
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repelled  by  the  positively  charged  grid.  Unless  otherwise  stated,  all  calculations  using  current  density  from  the  nude 
Faraday  probe  use  corrected  values  beyond  30  degrees  off  thruster  centerline.  Beam  divergence  calculations  with 
the  gridded  Faraday  probe  data  use  uncorrected  values. 


IV.  Experimental  Results  and  Discussion 

Thruster  efficiency  and  the  various  utilizations  are  analyzed  along  a  line  of  constant  power  (6  kW)  and  anode 
flow  rate  (20  mg/s)  to  study  the  losses  that  dominate  each  mode  of  operation.  Propellant  utilization  can  be 
calculated  using  the  experimentally  determined  thruster  efficiency  along  with  voltage  utilization  from  the  RPA  and 
current  utilization  from  the  Faraday  probe.  The  primary  drawback  to  this  approach  is  the  error  associated  with 
integrating  ion  current  density  in  the  plume  to  calculate  the  current  utilization.  An  artificially  large  beam  current 
would  decrease  propellant  utilization. 

The  anode  thruster  efficiency,  energy  efficiency,  and  three  utilization  efficiencies  are  shown  in  Fig.  8. 
Propellant  utilization  is  the  driving  loss  mechanism  for  both  the  constant  flow  and  constant  flow  rate  conditions.  A 
low  propellant  utilization  can  be  attributed  to  low  ionization  fraction  and/or  a  large  momentum  divergence  angle. 
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Figure  8.  Thruster  efficiencies  for  constant  anode  flow  rate  of  20  mg/s  (5a)  and  constant  power  of  6  kW  (5b) 
with  increasing  discharge  voltage. 

In  both  plots,  the  voltage  and  current  utilization 
are  nearly  equal  above  200  V  discharge.  At 
discharge  voltages  below  200  V,  the  potential 
required  for  ionization  becomes  a  significant 
fraction  of  the  anode  acceleration  potential.  This 
loss  of  available  acceleration  potential  is  seen  in  the 
reduced  voltage  utilization  in  the  20  mg/s  case  for 
100  Vand  150  V. 

An  inspection  of  the  relationship  between 
thruster  efficiency,  energy  efficiency,  and  propellant 
utilization  in  Fig.  9  reveals  similar  characteristics 
between  the  lines  of  constant  power  and  flow  rate. 

For  discharge  voltages  >300  V,  the  propellant 
utilization  and  energy  efficiency  are  nearly 
equivalent.  However,  at  400  V  discharge,  the  line 
of  constant  anode  flow  continues  the  previous  trend 
while  the  line  of  constant  power  deviates  to  lower 
energy  efficiency  with  higher  propellant  utilization. 

A  possible  mechanism  of  this  phenomenon  will  be 
revisited.  Since  discharge  voltage  appears  to  be  the 


Energy  Efficiency 


Figure  9.  Relationships  between  propellant  utilization, 
thruster  efficiency,  and  energy  efficiency  for  a  line  of 
constant  anode  flow  (20  mg/s)  and  constant  power  (6  kW). 
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primary  efficiency  driver,  further  analysis  will  focus  on  variations  with  applied  anode  potential. 

The  average  particle  charge  can  be  evaluated  using  the  quantity  (l-r)/E2  =  /  =  f  Q,  with  profiles  of  constant 
ionization  fraction.  Although  ionization  fraction  will  most  likely  vary  over  the  wide  range  of  discharge  voltages 
studied,  the  traces  in  Fig.  10  bound  the  likely  range  of  average  charge.  In  the  20  mg/s  case  shown  in  Fig.  6a,  the 
average  ion  charge  increases  exponentially  lip  to  200  V  discharge  and  plateaus  at  Q~  1 . 3  for  the  case  of  95% 
ionization.  While  it  is  possible  average  charge  varies  for  discharge  voltages  >200  V,  the  ionization  fraction  would 
vary  by  the  same  amount  since  the  /-Factor  is  constant.  The  lower  limit  on  Q  is  given  by  the  line  of  100% 
ionization,  where  f;=l  and  /=Q.  Thus,  for  discharge  voltages  >200  V,  the  average  charge  must  be  >1.23,  which 
implies  a  significant  fraction  of  multiply  charged  ions  in  the  thruster  discharge.  For  constant  power  operation  from 
200  -  400  V  in  Fig.  10b,  the  average  charge  decreases  with  increasing  voltage.  Similar  arguments  for  the  minimum 
allowable  Q  hold,  where  the  average  charge  must  be  greater  than  or  equal  to  the  line  of  constant  f;=l. 


10a) 


10b) 


Figure  10.  Average  ion  charge  for  constant  anode  flow  rate  of  20  mg/s  (10a)  and  constant  power  of  6  kW 
(10b)  with  increasing  discharge  voltage. 

The  maximum  momentum  divergence  half-angle  can  be  examined  using  the  values  of  the  Ei,  (1-P),  and  y  to 
create  similar  profiles  for  lines  of  constant  ionization  fraction,  as  shown  in  Eq.  (16)  and  (17).  For  an  isotropic 
distribution  of  divergences,  the  divergence  component  of  the  propellant  utilization  is  the  factor  by  which  thrust  of  an 
ideal  1 -dimensional  rocket  jet  is  reduced  to  its  2-dimensional  divergence. 
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Profiles  of  maximum  momentum  divergence  half-angle  are  shown  in  Fig.  11.  Calculated  divergence  of  the  20  mg/s 
conditions  shown  in  Fig.  11a  display  an  exponential  decrease  with  increasing  discharge  voltage,  and  exposes  beam 
divergence  as  a  leading  loss  mechanism  at  low  discharge  voltage.  In  this  case,  the  line  of  fi=l  places  an  upper  limit 
to  the  momentum  divergence  half-angle  of  ~35°  at  300  V.  In  Fig.  lib,  lines  of  constant  power  show  a  linear 
relationship  with  discharge  voltage,  where  the  maximum  momentum  divergence  angle  decreases  by  10°  from  200  - 
400  V  for  constant  6  kW  operation. 

A  comparison  between  the  charge  divergence  angle  and  maximum  momentum  divergence  angle  is  displayed  in 
Fig.  12.  The  90%  and  95%  beam  divergence  angles  calculated  from  Faraday  probe  sweeps  of  current  density  in  the 
far-field  plume  are  higher  than  the  momentum  divergence  angle  for  all  discharge  voltages,  with  the  lone  exception 
of  300  V,  20  mg/s.  The  most  likely  cause  of  this  discrepancy  is  the  presence  of  low  energy  CEX  ions  on  the  thruster 
wings  contributing  a  significant  fraction  to  the  integrated  charge  divergence.  It  is  expected  that  these  low  energy 
ions  have  little  effect  on  the  thrust,  and  may  supply  an  artificial  increase  to  the  calculated  beam  current  and  charge 
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Figure  11.  Momentum  divergence  half-angle  for  constant  anode  flow  rate  of  20  mg/s  (11a)  and  constant 
power  of  6  kW  (lib)  with  increasing  discharge  voltage. 
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Figure  12.  Comparison  of  maximum  momentum  divergence  half-angle  with  charge  divergence  half-angle 
for  lines  of  constant  anode  flow  rate  (20  mg/s)  and  constant  power  (6  kW)  with  variations  in  discharge  voltage 
(12a)  and  discharge  current  (12b). 

divergence  angles.  In  this  case,  the  overestimate  in  beam  current  would  lower  propellant  utilization  and  raise  the  %- 
Factor,  which  in  turn  would  cause  a  reduction  in  the  momentum  divergence  calculation  and  an  increase  in  the 
average  charge. 

Looking  at  the  limits  of  specific  terms  also  aids  in  further  bracketing  the  plasma  parameters.  For  the  point  at  400 
V  and  6  kW,  the  y- Factor  is  0.18  and  the  quantity  <bI)lvfi  =  0.88.  Using  these  values  and  the  limits  of  each  term 
(Q>1,  1  >  fi,  1  >  ®Div),  a  table  of  values  can  be  created  that  displays  probable  plasma  properties. 


Table  1:  Parameter  Sets  of  Plasma  Properties  for  400  V  discharge. 
400  V,  6  kW _  400  V,  20  mg/s 


ft 

Q 

<Ddiv 

a 

0.88 

1.34 

1.00 

0.00 

0.90 

1.31 

0.98 

12.14 

0.92 

1.28 

0.96 

17.05 

0.94 

1.26 

0.94 

20.75 

0.96 

1.23 

0.92 

23.81 

0.98 

1.20 

0.90 

26.46 

1.00 

1.18 

0.88 

28.82 

fl 

Q 

(Pdiv 

a 

0.88 

1.40 

0.95 

17.44 

0.90 

1.37 

0.93 

21.22 

0.92 

1.34 

0.91 

24.35 

0.94 

1.31 

0.89 

27.05 

0.96 

1.28 

0.88 

29.44 

0.98 

1.26 

0.86 

31.61 

1.00 

1.23 

0.84 

33.59 
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Highlighted  cells  show  the  most  probable  sets  of  consistent  plasma  properties  and  provide  further  insight  with  a 
minimal  number  of  diagnostics.  An  analogous  table  can  be  created  to  compare  the  change  created  by  increasing 
current  from  the  previous  example  of  400  V,  6  kW  to  the  400  V,  20  mg/s  conditions.  Analyzing  the  plausible 
plasma  properties,  it  appears  the  average  charge  Q  increases  with  discharge  current.  In  addition,  depending  on  the 
increase  in  Q,  the  momentum  divergence  angle  may  also  have  increased. 


V,  Summary  and  Future  Work 

A  systematic  study  of  Hall  thruster  performance  and  far-field  plume  measurements  from  100  to  400  V  discharge 
demonstrated  the  extraction  of  plasma  properties  with  an  analytical  efficiency  architecture.  This  method  accounts 
for  the  effects  of  the  ion  velocity  distribution  function  (VDF),  electron  recycle  fraction,  plume  divergence,  and  the 
fractional  loss  of  acceleration  voltage  within  the  discharge  chamber.  Propellant  utilization  was  the  driving  loss 
mechanism  for  both  the  constant  flow  and  constant  flow  rate  conditions.  A  low  propellant  utilization  can  be 
attributed  to  low  ionization  fraction  and/or  a  large  momentum  divergence  angle.  Physical  limits  of  plasma 
properties  enabled  limits  to  be  set  on  consistent  sets  of  performance  parameters,  including  the  average  ion  charge, 
ionization  fraction,  and  momentum  divergence  half-angle.  Including  ExB  measurements  or  LIF  measurements  with 
this  technique  would  close  the  set  of  equations  and  enable  the  calculation  of  a  full  set  of  discharge  parameters, 
including  the  aforementioned  properties.  Internal  and  near-field  plasma  potential  and  density  measurements  will  be 
conducted  to  further  understand  the  local  mechanisms  affecting  Hall  thruster  operation. 
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Appendix 


Table  1 :  Thruster  telemetry  data  for  lines  of  constant  mass  flow  (20  mg/s)  and  power  (6  kW). 


Test 

Point 

VD  (V) 

Id  (A) 

Pd  (A) 

Mass  Flow  (mg/s) 

Cathode 
Potential  (V) 

Background 
Pressure  (Torr) 

Thrust 

(mN) 

Anode 

Cathode 

8 

100 

22 

2245 

20.2 

1.4 

-7.4 

1.4E-05 

181 

9 

150 

22 

3340 

20.2 

1.4 

-9.1 

1.4E-05 

263 

10 

201 

21 

4289 

20.2 

1.4 

-10.2 

1.3E-05 

318 

11 

301 

20 

5997 

20.2 

1.4 

-11.6 

1.2E-05 

405 

13 

401 

20 

8040 

20.2 

1.4 

-10.7 

1.2E-05 

475 

19 

400 

15 

6044 

15.8 

1.1 

-13.3 

9.2E-06 

366 

20 

200 

30 

6031 

27.3 

1.9 

-10.3 

1.9E-05 

435 

21 

600 

10 

6029 

10.6 

0.7 

-18.1 

7.1E-06 

304 

Table  2:  Experimental  Parameter  Groups,  Utilization  Efficiencies,  and  Plasma  Parameters  (for  f;=l). 


Test 

Experimental 

(1-P) 

(1-r) 

<D 

Thrust 

For  f,  =  1 

Point 

E, 

E, 

Efficiency 

X 

OvDF 

Odiv 

a  (Deg) 

8 

0.54 

0.66 

0.80 

0.73 

0.62 

0.36 

1.10 

1.10 

0.99 

0.63 

54 

9 

0.76 

0.67 

0.87 

0.81 

0.73 

0.51 

1.21 

1.21 

0.99 

0.74 

44 

10 

0.84 

0.70 

0.87 

0.86 

0.78 

0.58 

1.24 

1.24 

0.99 

0.78 

40 

11 

0.91 

0.75 

0.90 

0.92 

0.82 

0.68 

1.23 

1.23 

0.99 

0.83 

35 

13 

0.94 

0.74 

0.92 

0.91 

0.83 

0.69 

1.23 

1.23 

0.99 

0.84 

34 

19 

0.91 

0.77 

0.90 

0.91 

0.86 

0.70 

1.18 

1.18 

0.99 

0.87 

30 

20 

0.86 

0.66 

0.85 

0.87 

0.77 

0.57 

1.31 

1.31 

0.99 

0.78 

40 

21 

0.93 

0.78 

- 

- 

- 

- 

- 

- 

- 

- 

- 

Table  3:  Far- field  Plasma  Plume  Measurements  at  1 -meter  on  thruster  centerline. 


Test 

Point 

Nude  Faraday  95%  Flalf-Angle 

Nude  Faraday  90%  Half-Angle 

Integrated  Beam 
Current  (A) 

Plasma 
Potential  (V) 

Most  Probable 
Ion  Voltage  (V) 

Avg.  (deg) 

±  Deviation  (deg) 

Avg.  (deg) 

±  Deviation  (deg) 

8 

72.9 

1.1 

63.2 

1.1 

16.3 

8 

88 

9 

57.1 

0.9 

47.4 

0.8 

18 

12 

142 

10 

50.1 

1.0 

44.4 

2.2 

18.4 

13 

188 

11 

35.8 

3.5 

29.9 

2.6 

18.3 

14 

284 

13 

37.0 

1.4 

30.7 

1.7 

18.3 

15 

382 

19 

37.0 

1.4 

30.7 

1.7 

13.7 

17 

376 

20 

60.6 

0.5 

50.6 

0.4 

26.3 

12 

182 

21 

40.7 

2.8 

33.5 

2.4 

10.1 

- 

- 
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